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We use the variational matrix-product ansatz to study elementary excitations in the S — | ladder 
with additional diagonal coupling, equivalent to a single S = | chain with alternating exchange and 
next-nearest neighbor interaction. In absence of alternation the elementary excitation consists of two 
free S = \ particles ("spinons") which are solitons in the dimer order. When the nearest-neighbor 
exchange alternates, the "spinons" are confined into one S = 1 excitation being a soliton in the 
generalized string order. Variational results are found to be in a qualitative agreement with the 
exact diagonalization data for 24 spins. We argue that such an approach gives a reasonably good 
description in a wide range of the model parameters. 
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I. INTRODUCTION 

Spin ladders continue to attract much attention 
as structures intermediate between one- and two- 
dimensional ones and possibly important ipr the under- 
standing of the high-T c superconductivity.!!! On the other 
hand, there exists a close relationship between "gener- 
alized" spin ladders (with an additional diagonal cou- 
pling), antiferromagnetic chains with frustrating next- 
nearest neighbor interaction, and the Haldane systems. 

In the present paper we study elementary excitations 
of the generalized 5 = \ spin ladder model (equivalent 
to a single zigzag spin chain with alternation and frustra- 
tion). The model is described by the Hamiltonian (see 
Fig. @) 

H = J2 S lin S 2 ,„ + (1 + 7) S b»S 2 ,„ + i 
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+ A]T(S 
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where A > determines the strength of the next-nearest 
neighbor interaction, and 7 corresponds to alternation 
of the nearest-neighbor exchange, whose strength is set 
to be unity. This model has rich behavior depending 
on the values of parameters A and 7, and itjL_phase dia- 
gram (see Fig. ^) is rather well understoocloa (at least 
for the half-plane A > 0; for negative A the situation 
is less cleam). Without loss of generality, we will as- 
sume that 7 <, 0, since there is an obvious symmetry 
transformationB relating the half-plane 7 > with the 
strip — 1 < 7 < 0. The "symmetry line" 7 = is pecu- 
liar because it is the line of transition between dimerized 
phases with different signs of the dimer order. Within 
the interval < A < A c w 0.24 this transition is of the 
second order, the ground state is unique and nondimer- 
ized, and the corresponding spectra are gapless; on the 
other part of this line the transition is of the first or- 
der, so that for A > A c , 7 = there are two degenerate 
dimerized ground states, and the system is gapped. The 



transition at A = A c is well studied, see Ref. I for a re- 
view. Everywhere except at the symmetry line the model 
(|l|) has unique ground state with a finite gap above it. 
In the limit 7-^—00 the diagonal spins form effective 
5 = 1 units, and the system becomes equivalent to the 
5=1 Haldane chain, with the effective coupling constant 
(l + 2A)/4. 

Elementary excitations of the generalized ladder, how- 
ever, have been studied to much less extent than its 
ground state properties. The Heisenberg point (A = 0, 
7 = 0) is exactly solvable by means of the Bethe ansatz 
technique, and the elementary excitations, are pairs of 
noninteracting S = | entities ( "spinons" ) The ground 
state contains the Fermi sea of spinons, and the excita- 
tions are of the particle-hole type. 

It is also known that at the so-called Majumdar- Ghosh 
(MG) point (7 = 0, A = h), where the exact (twofold 
degenerate) ground state is a simple product of singlet 
dimers,El the elementary excitation can be approximately 
constructed as a pair oLunbound spins above the com- 
pletely dimerized stateE The elementary excitation is 
composed of two S = h entities which are kinks in the 
dimer order and resemble spinons in that they are "al- 
most free" (i.e., form scattering states for most values of 
momenta) in case of unbroken translational invariance. 
When nonzero alternation is present, one may expect 
that those "spinons" are confined into a single 5 = 1 
particle. 

In the Haldane limit (infinite alternation) the system 
has long-range hidden kstring) order,u and the elemen- 
tary excitation is knownll3il3 to be a soliton in the string 
order. Ihexoncept of string order was generalized to spin 
Iadderal3il3 and it was shown that several ladder r»ed- 
els exhibit long-range generalized string correlations.^ It 
was arguedlj that the gapped phase of the spin ladder 
is the same as the Haldane phase of the effective 5 = 1 
chain. 

We construct the variational ansatz for the wave func- 
tion of the elementary excitation in the form of a matrix 
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product, using the recently proposedu matrix-product- 
states approach to the description of the ground state 
properties of spin ladders. The ground state in this ap- 
proach has a built-in generalized string order, and we 
construct the excitation as the 5=1 composite par- 
ticle being a kink in the string order and consisting of 
two bound S = \ entities. The wave function of the kink 
pair contains only one variational parameter £ having the 
meaning of the average pair size (the localization length) . 
We show that at the MG point £ is infinite and our ansatz 
reduces to that of Shastry and Sutherland;!! when mov- 
ing away from the MG point, £ drops down very quickly 
to the value of about a few lattice constants. At the 
"regular" ladder point (7 = — 1, A = 1) the dispersion 
relation as obtained from our ansatz agrees well with the 
results of other authors using different technics (see Ref. 
[j] and references therein). In the Haldane limit £ goes to 
zero, and our wave function transforms into -the "crack- 
ion" ansatz introduced by Fath and SolyorrJlil (see also 
Ref. |l6|) for the description of the Haldane triplet in the 
Affleck-Kennedy-Lieb-Tasaki (AKLT) model. We com- 
pare our variational results with the numerical data ob- 
tained through the exact diagonalization of a finite (24 
spins) ladder system, and find a reasonable agreement 
between the two approaches. We conclude that our sim- 
ple variational ansatz allows us to study analytically at a 
qualitative level the crossover from free to strongly bound 
spinons, giving a reasonably good description in a wide 
range of physical models. 

The paper is organized as follows: in Sect. |l|we intro- 
duce our ansatz for the elementary excitation. In Section 
III we present results from the variational calculation in 
comparison with numerical data. The elementary exci- 
tations in different regions of the phase diagram are dis- 
cussed. Finally, Sect. IV contains concluding remarks. 



II. TWO-SPINON "COMPOSITE PARTICLE" 
ANSATZ 

Recently,! a variational wave function for the descrip- 
tion of the ground state properties of generalized spin 
ladders was proposed in the form of a matrix-product 
(MP) state. ThcJVIP representation was first discussed 
by Fannes et al£3 in an abstract manner and later by 
Klumper et al$3 for the S = 1 deformed VBS chain, 
and has found since then numerous applications in exact 
and variational calculationsE3t3a. For periodic bound- 
ary conditions, the trial wave function for the ground 
state of the ladder consisting of 2N spins can be written 
as 

|*o) = Tr (g 1 g 2 ■ ■ ■ g N ) , where 

gi(u,v) =u(l- |s)ij +v ^2 ^> ' ( 2 ) 

Here the elementary matrix g, is constructed from the 
singlet state \s)i and three triplet states of the lad- 



der diagonals; are the Pauli matrices in a spherical 
basis, (To = <J Zl o±\ = =Fl/V2 (a x ± i<J y ). The parame- 
ters u and v in case of absence of external magnetic field 
can be chosen real and satisfy the normalization condi- 
tion u 2 + 3v 2 = 1. 

The wave function (|2|) has the following remarkable 
propertieslj (i) for arbitrary u and v it is a global sin- 
glet (see also Ref. ^2] for details); (ii) it has a built- 
in generalized string order defined as "diluted antiferro- 
magnetic order": and \t-i) should occur in perfect 

antifcrromagnetic sequence, arbitrarily diluted by |s)'s 
and |io)'s; (hi) both degenerate dimer ground states at 
the Majumdar- Ghosh point can be written in the above 
form, as well as the "valence-bond" AKLT stateEl which 
approximates the ground state in the effective S = 1 
Haldane limit. The two dimer ground states at the MG 
point correspond to u = 1, v = and u = v = \ respec- 
tively, and the AKLT state corresponds to zero singlet 
weight [u = 0, v = l/\/3). The state (|j) has very short- 
ranged correlations and therefore cannot be considered as 
a good approximation in the gapless region of the phase 
diagram; we will thus restrict ourselves to the study of 
the gapped phase only. 

The variational energy Eq — (^ol^l^o) calculated 
with the trial function (^|) has at most two minima as 
a function of u, v: one is always located at u = I, v = 
and corresponds to singlets on the diagonal links, and 
the position of the other minimum (ito, ^0) depends on A 
and 7; this latter minimum is absent in certain region of 
the phase diagram. The two minima have equal energies 
only at the MG point, and for any other choice of the 
model parameters they are inequivalenta 

We construct the trial wave function for the elemen- 
tary excitation, requiring the following: (i) it should be a 
triplet; (ii) it should be a soliton in the generalized string 
order as defined above; (iii) it should be able to repro- 
duce the ansatz of Shastry and Sutherland for the MG 
pointj3 i.e. a pair of unbound spins connecting two de- 
generate dimer ground states, and the "crackion" ansatz 
of Fath and Solyom.M One can check that the following 
construction satisfies the above requirements: 

n n N 

|n,n'; Ai )=Tr{(n^ ( ° ) )-i II « II ^} ■ (3) 

i— n+1 i—n' + l 



i=l 



Here g^ — <7i(uo, «o) denotes the matrix corresponding 
to the variational ground state, and cji = gi(u = 1, v = 0) 
is the matrix describingsinglet dimers on the diagonals 
(it is easy to see from (||) that gi is proportional to the 
unit matrix). In fact, the state (H) describes two domain 
walls which correspond to transitions between the two in- 
equivalent variational ground states mentioned above; it 
is worthwhile to remark that the states |n, n'; /x) are not 
mutually orthogonal. The presence of cr^ ensures that 
this state breaks the generalized string order and has the 
total spin S = 1 and its z-projection S z — fi (the general 
technique of constructing MP states with given quantum 
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numbers S, S z can be found in Ref. g2|). The structure 
of the ansatz (||) is schematically shown in Fig. ||. 

At the MG point uq = v q = h , and one can straightfor- 
wardly check that in this case Eq. (^) describes a pair of 
unbound spins ^ separating two completely dimerized re- 
gions, i.e. it reduces to the ansatz of Shastry and Suther- 
land. On the other hand, if n = n' and uq = 0, then the 
state (||) is exactly the same as the "crackion" of Fath 
and Solyom. 

We use the following (unnormalized) trial wave func- 
tion of the composite two-spinon excitation with a given 
total momentum k: 



k(n+n')/2 e iq(n-n')/2 e -(n'-n)/t\ n ^ n > 



(4) 



It contains two variational parameters q and £ which can 
be considered as real and imaginary part of the relative 
momentum of the two spin-i entities forming our com- 
posite particle. The parameter £ has the meaning of a 
localization length, or the average size of the compos- 
ite object, and nonzero q corresponds to the excitation 
of some internal degree of freedom. If the localization 
length diverges, £ — > oo, the wave function describes a 
triplet scattering state of two spinons, and finite £ corre- 
sponds to a bound state. 

On the disorder line 7 = 2 A — 1, according to Ref. |^, 
the parameters Uq and v both equal |, and the struc- 
ture of our variational ansatz becomes rather obvious: it 
describes a bound state of two Shastry-Suthcrland kinks 
("free spins" in a completely dimerized chain) with the 
localization length £. 

The energy for such an excitation can be calculated in 
the usual way: 



(k,q;n\H - E \k,q;^) 
(k;n\k;n) 



(•5) 



and has to be minimized over £ and q, separately for any 
given k (as we will see below, it turns out that optimal 
£ strongly depends on k). In this way one looks for a 
lowest variational state in a subspace with the total spin 
5=1 and certain momentum k. Calculating the av- 
erages in (|^) involving MP states_can be done with the 
help of the standard technique.EaO The final expression 
for E(q,£,k) is quite cumbersome because of the com- 
plicated structure of our trial wave function (|^), so that 
we present here only the resulting dispersion plots for a 
number of representative points of the phase diagram. 
The minimization has been performed numerically. 

Another, simpler ansatz can be obtained if one forces 
([|) to be strongly localized, i.e. £ — > 0. Then only the 
configurations with n = n' survive, and we obtain 

\k; (j) = e lkn \n; fj,} , where 

n 

k M) = Tr {(ffW • • • 9^Hi&& -gf}. (6) 



Such an ansatz may be called "regular," or localized 
crackion (in contrast to our two-spinon ansatz (^) de- 
scribing an "extended crackion" ) because it exactly re- 
produces the structure of the wave function proposed by 
Fath and Solyom, and the only difference is in the gener- 
alized concept of the string order (or, in other words, in 

the fact that matrices gf 1 are allowed to contain sin- 
glet states). This "regular" crackion ansatz is essen- 
tially equivalent tq-the construction introduced recently 
by Nakamura et alt3 within a slightly different approach 
using the Kennedy- Tasaki unitary transformation. 



III. VARIATIONAL RESULTS FOR THE 
EXCITATIONS AND COMPARISON WITH 
NUMERICAL DATA 



In this section we present variational and numerical re- 
sults for dispersion relations, lowest modes and their cor- 
responding wave vectors for various points of the phase 
diagram. In the first subsection we study the spinon-type 
excitations on the symmetry line 7 = in the vicinity of 
the MG point. Variational energies for scattering states 
and bound states are computed at the MG point A = i, 
and the difference in behavior of the spectrum at A > | 
and A < i is discussed. In the second subsection we 
investigate the consequences of alternation. When one 
moves off the symmetry line, the crossover from spinon- 
type excitations to ladder- type ("crackion") excitations 
occurs, being characterized by the change of the localiza- 
tion length £(fco) from infinity to zero; here ko denotes 
the wave vector of the lowest energy mode. We show 
that this crossover takes place at the ko — tt boundary of 
the incommensurate region where ko changes gradually 
from to 7r. These results are compared with numerical 
data from exact diagonalization for 24 spins. Because 
the correlation lengths for the points taken into account 
are rather small, the numerical data are very close to the 
thermodynamic limit and therefore we do not perform 
any finite-size extrapolation. 



A. Elementary excitations on the symmetry line 
7 = 



The ansatz (|3j,[|) obviously becomes inadequate close to 
the symmetry line 7 = (except for the MG point), be- 
cause two variational minima of the ground state energy 
are inequivalent while the true ground state is twofold de- 
generate on this line. However, the interval of the sym- 
metry line in the vicinity of the MG point (i.e. 7 = 
and A close to |) can be studied with the help of the 
Shastry-Sutherland-type ansatz: it is sufficient to put 
u = v = 1/2 in (||,||). In the limit £ — * 00 one gets 
exactly the ansatz of Ref. || and can calculate variational 
energies for the scattering states. On the other hand, if 
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we do not force £ = oo, the energy of the lowest bound 
state can be calculated for each value of the total momen- 
tum fc. For the scattering states it is possible to obtain 
a compact analytical expression for the energy of such a 
two-particle excitation: 



E(M) = £ (^) +£ (i i i) 



e(fc) = -(5 + 4cosfc) + (1 



(7) 




5 cos k + 4 
5 + 4 cos k 



here e(k) has the meaning of the spinon dispersion. At 
A = | it coincides with the result of Ref. ||. (It should be 
mentioned that since in this paper we consider the gen- 
eral case of an alternated chain, momenta in Eqs. (|4|]7]) 
are defined in the halved Brillouin zone, in contrast to 
Ref. |^). The dispersion relation is determined by the 
lower boundary of the two-particle continuum. The re- 
sulting A-dependence of the gap A is shown in Fig. ||, tp» 
gether with the numerical results by White and Affieck.E3 
One can see that quantitatively this approach yields rea- 
sonable results only in the close vicinity of the MG point; 
nevertheless, at the qualitative level it correctly predicts 
closing of the gap when decreasing A and the existence 
of a maximum in -Eg (A) at A slightly greater than i. 

Characteristic plots of the dispersion of scattering 
states in the vicinity of the MG point are shown in Fig. 
||. From (0) one can see that for A less than = ^ the 
single-spinon dispersion has a minimum at k = fco = n, 
and for larger A this minimum shifts towards fco lower 
than 7T. Thus, the lower boundary of the two-spinon con- 
tinuum always has the minimum at fc = (corresponding 
to q = 2fco), and for A > A^ there appears another mini- 
mum at fc = fco = 2(tt — fco) (which corresponds to q = 0). 
When A increases, this second minimum gets more pro- 
nounced. 

Appearance of the lowest mode with incommensurate 
wave vector is closely xelated to the existence of the so- 
called disorder pointsEjEE where spin-spin correlations 
in real space become incommensurate. Strictly speaking, 
the point where the wave vector of the lowest mode be- 
comes incommensurate corresponds not to the disorder 
point itself, but to the so-called Lifshitz point where the 
correlation function peak in momentum space (i.e., peak 
in the structure factor S(q)) starts moving from com- 
mensurate to incommensurate q. Generally, the Lifshitz 
point does not coincide with the disorder point and is 
situated at some small distance from the boundary of 
the incommensurate region. In our variational calcula- 
tion the wave vector starts to change at A w = while 
the disordcx,linc (line of disorder points) is numerically 
establishedoH to be 7 = 2A — 1, which means that at 
7 = the disorder point is A = h < , in agreement 
with the above. 

We would like to end this subsection by pointing out 
the role of bound states. At the MG point it is knowntj 
that bound states are lower in energy than scattering 
states for wave vectors fc close to the zone boundary. This 



was obtained in Ref. || by using the variational estimate 
for the upper bound of the dispersion (see below). We can 
capture the dispersion of the lowest bound state in our 
approach if we minimize with respect to £ for each fc. For 
values of the total momentum 0.687T < fc < n we obtain 
> as shown in Fig. The same feature, namely 
the appearance of bound states as the lowest energy ex- 
citations, can be expected for A > A^ and fc around the 
midpoint in between the two dispersion minima, but we 
did not perform this calculation. Generally, more than 
one bound state may exist p3 unfortunately, within the 
present approach one can access only the lowest bound 
state. 



B. Dispersion relations and crossover from loosely 
bound to tightly bound spinons 

First of all, we would like to focus on the disorder line 
7 = 2 A — 1. As soon as one moves off the MG point 
(A = i, 7 = 0) towards the dimer point (A = 0, 7 = — 1), 
the energy of the MG ground state with singlets on diag- 
onals picks up an energy proportional to the size of the 
system, while the other MG ground state with singlets 
on the rungs remains the true ground state. As a conse- 
quence, only bound states of spinons can survive: this is 
a typical confinement situation. The wave vector of the 
lowest mode is still fc = 0. 

The main feature of this excitation is the fc-dependent 
derealization: The ground state is formed by a prod- 
uct of rung singlets. If we now replace one singlet by a 
triplet and superpose with wave vector 7r, we obtain an 
exact eigenstate on the disorder line- whose energy gives 
the upper bound of the dispersion EH To obtain the low- 
est mode we delocalize the up spins, which make up the 
triplet, with the amplitude e~^l r as shown in Fig. ^ and 
superpose with the wave vector fc = 0. The bandwidth 
for this excitation is illustrated in Fig. |7]; one can see 
that the bandwidth increases along with the increase of 
the localization length £ from to 00 on the way from the 
dimer point to the MG point (here £ is not to be confused 
with the spin correlation length!). At the MG point we 
reproduce the result of Shastry and Sutherland for the 
gap A = 1/4. One may say that here we observe the 
crossover from loosely bound spinons ("extended crack- 
ion", 1 <C £ < 00) to spinons tightly bound into the 
Haldane triplet ("crackion," £ <C 1), even though the 
crackion at the dimer point is a trivial (dispersionless) 
excitation. 

The same physical picture of crossover should be valid 
for any path beginning somewhere at the symmetry line 
and ending somewhere at sufficiently negative 7; on the 
symmetry line £(fco) should be infinite (fco denotes the 
lowest mode wave vector), and it decreases to zero when 
—7 is large enough. In the present approach we can 
observe this £ — > 00 behavior only for the MG point, 
because for any other point on the symmetry line our 
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two variational ground states become inequivalent and 
we lose the feature of twofold degeneracy. However, we 
believe that our approach remains reasonable for points 
which are far enough from the symmetry line, where the 
translational symmetry is explicitly broken and this ef- 
fect overrides the built-in dimerization of our variational 
ansatz. 

It is worthwhile to make a few remarks concerning the 
behavior of the real part of relative momentum q. Be- 
cause of the confinement, the energy of states with inter- 
nal motion, i.e. with g ^ 0, 2tt, has to be much higher 
compared to the states with zero relative momentum. 
In our variational calculations we were able to detect 
only one minimum which always occured at q = or 
2tt (strictly speaking, q = for tt < k < 2tt and q = 2tt 
for < k < tt, so that sets with q = and q = 2tt 
correspond to physically equivalent states). 

The crossover from 'extended crackion' to crackion is 
related to the change of wave vector ko of the lowest 
mode, due to the following remarkable feature which we 
observed from our calculations and which is illustrated by 
Fig. |^: everywhere in the gapped region of the phase dia- 
gram the property lim^—nr £(fc) = holds. Thus, once ko 
has changed from to tt, we know that the lowest mode is 
a 'usual' crackion because £(tt) — 0, so that the crossover 
takes place at the ko = tt boundary of the incommen- 
surate region. In order to determine the boundaries of 
the incommensurate region, we used our variational ap- 
proach and compared the results with exact diagonal- 
ization data. In Fig. ||(a) dispersion relations for a few 
points on the vertical line A = ^ are presented. Numeri- 
cal data and variational results are in a good agreement 
even though the lowest wave vectors are at slightly dif- 
ferent positions (one should keep in mind that numerical 
dispersions for finite chains consist of a finite number 
of points). Fig. ||(b) illustrates the change of ko when 
crossing the disorder line. Similarly to the situation on 
the symmetry line as described in the previous subsec- 
tion, one can see that ko starts to change from not 
exactly at the disorder line but slightly above it (i.e., the 
Lifshitz line is not identical to the disorder line). The 
comparison with the numerical data as is shown in Fig. 
^(a) confirms this property, within the numerical accu- 
racy (exact diagonalization of 24 spins leads to 12 values 
for the wave vector which is not sufficient to mark the 
incommensurate region precisely but allows qualitative 
comparison) . The boundaries of the incommensurate re- 
gion obtained by the variational calculation are presented 
in Fig. ^. It should be pointed out that our result for the 
fcor-= tt boundary does not agree with that of Pati et 
al.Ea For example, we obtain that the ko = tt boundary 
goes through the dimer point 7 = — 1, A = 0, while the 
corresponding curve C of Fig. 2 from Ref. |3^ crosses the 
7 = — I line at A ss 0.6. At present we cannot comment 
on the origin of this strong discrepancy. 

Finally, we would like to discuss the vicinity of the line 
7 = — I which includes the experimentally relevant ladder 
point (A = 1, 7 = —I). Figure O shows the dispersion 



curves for the "regular" crackion, the "extended" crack- 
ion and the numerical data for two points on this line 
(one should mention that our data agree rather well with 
those of Ref. 33). The dispersion curve of the extended 
crackion is located slightly below the curve of the crack- 
ion and coincides with it for k = tt, in agreement with the 
general property Unift^ £(/c) = (see also Fig. ||). One 
can also observe that the k = gap is slightly larger than 
2E(k = 7r), which indicates the repulsive character of the 
effective interaction between the elementary excitations 
in the ladder. Fig. [ll] displays the lowest crackion mode 
(k = tt) in comparison with exact diagonalization data 
along the horizontal line 7 = —1, A = 0...1. We con- 
clude that the ladder excitations can be described by the 
usual crackion ansatz, and the isotropic spin ladder has 
the same type of elementary excitation as the effective 
5 = 1 chain which appears as the limit 7 — > — 00 of the 
generalized frustrated ladder. In that limit the localiza- 
tion length £ collapses for all values of k because creation 
of singlets on the diagonals would cost infinite energy. 



IV. CONCLUSION 

We have presented a variational matrix-product ansatz 
for elementary excitations in the gapped phase of the 
S — i ladder with an additional frustrating diagonal cou- 
pling 1 + 7, 7 < 0; the strength of interaction along the 
legs is A, and the interaction along the rungs is chosen to 
be unity. This system is equivalent to the antiferromag- 
netic spin-i zigzag chain with alternating exchange (the 
magnitude of alternation is proportional to 7) and next- 
nearest neighbor interaction A. Our ansatz describes a 
triplet state of two S = h entities ("spinons") and al- 
lows one to interpolate between free and bound spinons 
by varying the parameter £ which has the meaning of 
a localization length (average distance between spinons 
in the pair) . This state is, constructed to be a soliton 
in generalized string orderfl and in the limit 7 — * — 00 
of effective S = 1 Haldane chain it coincides with the 
"crackion" ansatz proposed by Fath and Solyomliil if the 
localization length £ — > 0; for that reason we call our 
ansatz an "extended crackion." The limit £ — > 00 leads 
to the two-particle excitation of Shastry and Sutherland,El 
which corresponds to free spinons existing in absence of 
the alternation (i.e., on the symmetry line 7 = 0). 

Using our variational ansatz, we calculated dispersion 
relations for various points in the phase diagram. These 
results were compared to exact numerical diagonalization 
data for 24 spins and showed a reasonable agreement of 
the two approaches. The variational parameter £ was de- 
termined separately for each value of the total momen- 
tum k; it turns out that £(fc) has nontrivial behavior, 
particularly the property limfc-^ £ = was numerically 
observed in the entire range of studied model parameters. 

We determined the boundaries of the incommensu- 
rate region (strictly speaking, the corresponding Lifshitz 
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lines) by locating the wave vector ko of the lowest mode: 
for 7 > 2 A — 1 the wave vector is pinned at kg = 0, 
slightly after crossing this line it starts to change from 
to 7r, and, finally, at some other line it gets again pinned 
at ko = ir (in terms of the full Brillouin zone of the chain 
this corresponds to the change from ir to ir/2). 

We show that in the interval between the symmetry 
line 7 = and the ko — n boundary of the incommensu- 
rate region the dispersion of elementary excitations is well 
described by our bound spinons ansatz. The crossover of 
the lowest mode from "extended crackion" to "localized" 
crackion (i.e., from finite £(fco) to £(&o) = 0) occurs at the 
k = 7r boundary of the incommensurate region, due to 
the above-mentioned property of the function £(fc). For 
the isotropic ladder point (7 = — 1, A = 1) a localized 
crackion ansatz with £ = is sufficient to describe the 
excitations. 

To conclude, we propose a simple ansatz providing a 
reasonably good description of the elementary excitations 
in the gapped phase of the frustrated S = ^ chain with 
alternation for a wide range of the model parameters. 
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FIG. 2. Phase diagram of the generalized ladder. The interval of the line 7 = for < A < A c ~ 0.2411 is gapless, and the 
other part of the diagram is gapped, k = and k — n boundaries of the incommensurate region are variational estimates for 
the corresponding Lifshitz lines determined from the variation of the wave vector k of the lowest excitation-4«ode. The disorder 
line, where spin-spin correlations in the real space become incommensurate, for k — boundary is knownt3i3 to be 7 = 2A — 1; 
for k — 7r boundary the exact position of the disorder line is unknown. 
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FIG. 3. The structure of the variational two-spinon ansatz Ovals show the location of matrices; solid ovals denote g^ ' 

and dashed ones denote g. Thick solid links represent singlet bonds (bonds inside dashed ovals are always pure dimer bonds, 
and bonds between the solid ovals become purely singlet only on the disorder line 7 = 2 A — 1). 




8 




0.0 



scattering states 



bound states 




71/2 

total momentum k 



FIG. 6. Variational result for the dispersion at the MG point 7 = 0, A = ~ (solid line) in comparison with the numerical 
data for 24 spins (diamonds); the insert shows the momentum dependence of the localization length £(fc) for the bound states. 
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FIG. 7. The bandwidth E(0) — E(tt) and the localization length £ of the lowest energy mode on the disorder line 7 = 2A — 1. 
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FIG. 8. Momentum dependence of the localization length £ for two points on the line 7 = — 1. Note that £(71-) = 0. 
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FIG. 9. Dispersion curves from variational ansatz in comparison with the numerical data from exact diagonalization for 24 
spins: (a) three points on the line A = ^; (b) three points on the line 7 = — 2A in the vicinity of its crossing with the disorder 
line 7 = 2A — 1 (numerical data are taken from Ref. H) . 
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FIG. 10. The dispersion of elementary excitations on the line 7 = — 1: (a) the ladder point A = 1; (b) A = 5. Variational 
results from "extended crackion" (£ = £ m in(fc)) and "crackion" (£=0) approaches are shown along with the numerical data 
(O) from exact diagonalization. The dash-dotted curve shows the lower boundary of the two-particle continuum for 'extended 
crackion' ansatz. 
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FIG. 11. The energy of the lowest mode E(k = n) along the line 7 = —1; solid line represents the variational result and 
diamonds (O) are numerical points from exact diagonalization for 24 spins. 
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